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This paper, an expansion and discussion of a eeries of lectures given 


by Dr. F. N. Scheubel cf the Technical Institute, Carmetact, Germany, was 
undertaker by the writer as a theeie project while working for the degree 
of Mester of Science ir Aeroneutical Fngineering at the University of Mich- 
igan, Arn Arbor, Michigan. These lectures were given at the University of 
Michigan during the fall of 1951 while Dr. Scheubel was visiting the Uni- 
versity at the invitetion of Dr. F. W. Conlon of the Acronautical Engineer- 
ing Departnent. 

Dr. Scheubel presented the material in six lectures, the firet three 
devoted to quesi-static stebilitv corceptse, the last three to their use in 
the solution of the cquations of motion. Hise linited time orever.ted a de- 
talled accounting of the assumptions irvolved as well as minute exrlenationse 
of the approximations made. It is the purvose here to verify ard expand on 
his presentation to a degree consictent with the limited scope of such a 
paper. 

Briefly, his materiel covered the following. His lectures were re- 
etricted to symmetric or longitudinal motion. Fe developed the longitudi- 
mal equatione of motion and the quasi-static stability criteria at equili- 
brium and at constant speed. Hs sclved the equations useing these quasi- 
stetic stability concezpte for both the pruzoid and short period modes for 
the stick-fixed case. An amplitude-phase relation for the two variables 
aoplicable to the modes was discussed. He then introduced the degree of 
freedom about the elevator hinge line and solved for the etick-free case. 
Finally the effects of an elevator impulse were discueced. 

This approach to dynamic longitu.inal mction differs from standard 
methods mainly in the quasi-static stedility concey;te ae regards their in- 
sertion in the sclutions tc the ecustiorns of motion. In this respect, the 
method may, es illustrated, be safely applied only to conventional air- 
craft, i.e., a rigid body, where the number of degrees cf freedom and thus 
the complexity of soluticn is restricted. fhe enalysis of the mction as- 
euming a phygoii and a ehort veriod oscilletion servee only to illustrate 
the hendling of the quasi-static concerts. The limitation of the method 


in respect to acceptable separetion of the motion into these modes wes 








The equ:tione of motior. es eetabdliched by the standerd approech of 


reference 1 are relegeted to the flight path direction or wind axee in 
order that the further discussion will bce based on common esesumptions 
end epproximations. Dr. Scheubel’s equations as devcloped for a svetem 
reéeletive to the flight path will then »e ccroared term by term with the 
etandard equations. From this point orward, i.¢., »bécirming with non- 
dimensionalizing the equatione, the presentetion is as Dr. Scheuvel de- 


veloved it. 








2, M 


0,6 ,G or 6, As ,Aw 


NOMENCLATURE 


mass of the aircraft 


initial velocities along x and z body 
axis 


perturbation velocities along x and z 
body axis 


angular velocity about y axis (body or 
wind ) 


air forces and moments . 


perturbation of attitude angle, flight 
path angle, and angle of attack 


initial attitude angle 


aircraft radius of gyration about y 
axis 


moment of inertia about y axis 
initial flight path velocity 
velocity perturbation along flight path 


velocity along flight path following per- 
turbation 


lift, drag and air moment 

aircraft weight 

initial thrust and drag 

elevator deflection 

lift coefficient for aircraft 

air density 

wing area 

fimcrart drag cotiiicreac 
non-dimensional velocity perturbation, 


ratio of perturbation to total flight 
path velocity 
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ac" 


adopted standard time wit 


mass density (See Appendix A) 


non-dimensional time 

initial moment about y axis 

tail length, distance from aircraft 

c.g. to center of pressure of horizon- 
tal tail 

conjugate complex root of characteristic 


equation. kK is real part, I is imaginary 
part 


constant coefficients of equations 
position of aircraft center of gravity 
from most forward part of aircraft in 
percent of wing chords 

position of neutral point of aircraft 
from most forward part of aircraft in 
percent of wind chords 

angle of attack of hori. ontal tail 
lift coefficient of horizontal tail 
area of horizontal tail 


angular frequency of mode of motion 


periods of phygoid (p) and rotary (r) 
modes of motion 


wave length of phygoid 


distance between the hinge line and cen- 
ter of gravity of the elevator 


mase of the elevator 


radiue of gyration of the elevator 


hinge moment coefficient of elevator 


mean aerodynamic chord of the horizon- 
tal tail 





us elevator impulse <2 
H, constant altitude reference 
h altitude change 


Subscript ~ indicates steady state values in the response discussion. 
Any other symbols used are either obvious or are locally defined for ease 


in handling equations. 





II. Development of the Method. 
A. The equations of longitudinal motion. 
The standard development in a perfectly general way es accomplish- 


ed in reference 1 result in the equations of longitudinal motion of the forn, 
A.l = m(a' + Wa) = Xu + xv + Xa - mg cos 11,0 ; ) 


A.2 mw - U,a) = Zu + Cs 28 - mg cin re49 


A.3 mio g-I @-M ul + Mw + Ma . 
These equations are relative to body axes and contain the follow- 
ing eeSumptions and approximatione: 

(a) Initial symmetric steady motion is assumed. 

(bd) The air reactions do not depend on the rates of change of the 
variables, U, V and W, or their integrals. 

(c) Second order and higher terms of the air reactions are neglected, 
i.e., only infinitesimal disturbances from an initial eteady motion are 
treated. 

(d) The aircraft has a plane of symmetey and the steady motion about 
which the disturbances occur is symmetrical with regard to that plane. 

(e) The disturbance initially imposed on the system is unenforced and 
the controls are locked. This rules out an initial couple, Mo and the vari- 
ation of M with $. 

The variables are seen to be u, wand q. These are illustrated in Figure 
l. The relation between the body axes and the wind axes is shown in Figure 

n 
ec. The body axes are primed. Angles are measured positive counterclockwise 


from the horizontal reference for }* and *@., from the wind line for ™. 


|? 
In regard to these figures it is seen that, 


(a) Since q is a velocity and 6 a displacement, 6 = fa dt = AQKO+APL 
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(bd) V, 1s the steady estate velocity and, ve = ee 


Then there is the definition, V = os + AV 


dv _d(Av) 
also, ae o ae 


(c) In essence u is now WV, and w does not exist. 
The inertia terms of the equations of motion parallel and perpendicular 
to the wind direction and about the center of gravity become, neglecting the 


mercia force due to linear acceleration pempendicular to the wind direcvic, 


d(V) 
at 
@) Yr gine 
for A gen m \ at 
) 
ra HK A 
for A.3: I a(aA + ¥) 
wy até 


The weight componer.ts relative to the wind axes become, 
along x: mz sin ¥ 
along z: “mg COS of 
and these are only effected by A®& - 
There remains qnly the air reaction derivatives. The force created 
on the tail due to q ie the largest resulting from this disturbance and 
from experience it is known to be emall along the wind axes so and Z 


q 
eae neglected. 





In the standerd method the left and dreg contributed to the air reac- 
tion forces both elong the chord and normal to it. Aleo velocity pertur- 
betions existed along both axes. Using wind axes, cnly the velocity per- 
turbetion, AV, now exists. Any sinking velocity along the z wind axis is 
merely a change in angle of attack. The air reactions alcng x consist 
cnly of the drag in its entirety, neglecting variations of thrust. The 
dag varies with both & end V. The air reactions along z consist only’ 
of the lift in ite entirety and is affected by the same quantities. A | 
pure rotary perturbation about the center of gravity changes bcth & and 

* but only QO affects the linear forces. The same applies to the air mo- 
ments. The variations with &% eee ound by wind tunnel tests end any rota- 
tion of the model about the center of gravity ie considered pure % . This 
does not hold for the case of the effect on M of both A & and ao. These 
two perturbations make up @, which, when multiplied by the tail length, 
gives the effective sinking velocity of the horizontal tail. Thie sinking 
merocity in turn ic felt as a change in angle of attack of the tail. 


Thus the externel forces for A.1 become, 


oP 10 + 2 uy + SS (mg vin) at 


amd equetion A.l1 is, 
dite) “oD one 
nT Sq OK + Hp AV + mg cos FAY . 


hquétion A.2 is, 


ay OL OL . 
mV 7 7 5q 4% + 57 OV + me coin d AY 


rquation A.3 becomes, 


2. 
a(x + &)_ OM OM OM... . OM -; 
= Ly& / - : 
“yy age aa oe * By AV * Sa -% * BROT 


The use of the wind axes then hes the advantage of the aircraft's for- 


ward motion being along the x axis so that the lift and dreg forces always 








lie along these axes and resolution into normal end chord forces is not 
necegésary. Furthermore, no resolution of the velocity is required to be 
made alongs the axis of perpendicularity which is the case if the forces 
are resolved along the body axes. The disadvantage is thet these exes 
change contiruously sco that the momente and products of inertia change. 
For smell disturbances from initiel horicontal flight these discrepancies 
are considered negligible. Also the thrust does not seh) Quieas the wind 
emes for all aircraft and, for an indiviaiual aircraft, at all times. 

The continuity of this payer end the advantage of physical "feel" 
that is had from the development of the equations of motion from initial 
Suppoeiticn of wind axes can be best maintained by including euch a de- 


velopment in ite entirety at this point. This is the method used by Dr. 


Scheubel and appears to be in general use in Germany. See Keference 2. 





Figure 3. 


With the aircraft in an attitude end with the external forces eas 
Shown in Figure 3, the external forcee are equited to the time rate of change 
of momentum in the direction of the wind axis, i.e., along V, and one equa- 


ticn of equilibrium follows. 
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The externel forces are caused to chanse by emall changes or pertur- 
bations in the values of & , V, & and & . The forces in equilibrium 


then, following a small disturbance, are the initial values plus the in- 
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cremente due to the changes, ‘ 
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This is the result of a Taylor Series expension of the quantities T, t 
and W as functicne of A ,v, * and § in which second order terms and a- 
pove are neglected under the assumption that the changes ere smell. Assuming 
an initial condition of equilibrium, it is seen that this equation is compar- 
able to equation A.1 with one exception. Thue all the assumptions previously 
enumereted hold. The exception is the thrust term. It is easily seer that 
thrust is not affected by changes in KX , & , or §& . The etick-fixed case 
is considered here elso. Thus all terms note $ are neglected. 

This first equation is now manipulated into e form necessary for subd-, 
sequent solution. Initfel equilibrium and eseuming & smell gives, 


fo dD, - W sin eo = O Vande cos.4.= 1 or We EE 


e) v 
; 2 : 2 p 
also, L=C sfvs and, D=C,sfvs a 
| My “0 V there is 
Multiplying terms in T dv a end y and terms in D dy and + 
QO © 
W av Vor a VoD w oD 
- Ye aoe = . — — - : a he 
gz dt 7oT ov V °o D. Sv V oat — 
Ow, voor _8int. v aD. in D and letting a , 
tT, OV Gany 7D Ov ” Cuma om 
there is, 
Wav W aly. + AV) _w é(ay) 
edt ¢ At re at 
ee. We o (<2 ra) ln @ 7 0 in > ae Pat: os Ge Onn | y- ) Pe S 
DP > 8's" san d inv Sh we oh oad 2 
Oc 2 
— Tr : : Via -_ 
Using, the etanderc rotation, <aae “Dey end dividing bv I = S, Unerem. 
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Thie equation is non-dimeneioneli.ed in the ueual merner. The right 
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eo" a, ees ‘Sind ~ 5, Soay’ ° 
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cide of the ecuetion is non-dimersionel as it stends. Considering the left 


side, multiplying nunerator by V/V and the denominetor by t/t. 
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mecime unit, t_, is adopted by Tetting, z = = 1 then, 
é 6 Prves 8 


eel? oe 


y 
t = = = Cy - where the mass number, ie given by 
he fs g Y 


E pve l? ap 


Wy2/? iL ; 
LX > > - This maes number, or density factor, is 
spsve ap | wi 
discussed in Appendix A. Tow the non-dimensional time is given by, t/*. =¢ 
iG o~ . iy. 
60 that, d(s—) = dé and the varieble oV is transformed as, <a 
ine fir@®e equation of moticn eien becomesc, 
eee ee ee ee 
- Ge 6 ge ‘oe “o ‘6 Inv “bd, Sd inv 


feein referring to Figure 2, the emternal forces perpendicular to the 


wemd line aré equated to the centrifugal Borce or mas@ times the centriiuses 


P A "7 
mecelerstion. ul Vv 2 - V as L-Wcos & since, ul ee C, Deeevee 
QB ab g at | gz «at 


V == reeults in no chanze in velocity permencicular to the flight path and the 


linear acceleration in this direction is neglected. If a sm-1l disturbance 
is introduced, and, since lift does not vary with & «and the stick-fixed cone 
ig ascured, 


i aL aL 
rT —-— oo L ———— ye 4 sj ¢ aa ° 
ie Lo W cos ee = oO sy - W sin s tt 


oe (=; 


With the same acsumption as before, that of initial equilibrium, this 
equation isc ecen t> be the same as equ2ztion 4.2 and so the same assumptions 


and approximéetiors are implied. Fxpanding the derivatives, 
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= 5 RMI oo . ty r) 
au 2° eee * One 
vo 
aleo, W sin ws Sc, sir se so that the entire eau tion be- 
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YA dP, sou, (2) gest) uw a2) anos 
s Py St. tn cae Oo le ve a) 
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The left side is ron-dimensionsalized as before. Diviling the demom.nz tor 
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Dy t/t, tkeere ic, 
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So finallv there is the second equatim of action, 


: 3d ln cy 
= i Se a . : = 
ne x “TL (25; ae ET, U sie i, & “t ein {7 of or 


The final equation of the longitudinal motion is obtained by equeting the 


monerts to the momert of momentum about the center of grevity cf the aircraft. 


-~ 


eonsicer.ng Fleure 3 again, thie is, 


BA as at 
at : 


“4 


a 
So YY 
W .2 , 2 : 
where, a ey Is the moment of inertia end = ie the radius of «yration. ‘Since 
M ic a function of lift and dregs which in twrn are not affected vy & , nezlect- 
ing the influence of V, end with the etick-fixed esecum tior, there ie, 


a rS | 


Bee 2 (S py OU Ao 
> a F 
Ss 7 dt ye ae 
; OM OM OM . OM . 
a au 4 Sy - Soe ae” 


This is seen to be tre eame ae eqution “/3 essuming initial é@quilibriun. 
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Putting the change in moment in coefficient form, 





dc Oe ee 
M se 2, Mp eey2 ce 
cata a ros 7 oe pers ee ae PS at 
at 


. ~ a(aX ~ YY) ae 
since, AN +AX = at ee: 


An imoortant point to note here is the German vee of the equare root of 
the wing erea as the Baait veneel length term when converting moments to coef- 
ficient form. The reascning is that this quantity is more easily defineble 
for the variegated types of wings now in existence whereas the men aerodynamic 
chord is somewhat nebulous in definition in the litereture. Ancther reason 
ié the fact that the quantity r/s¥® ranges from C.8 to 1 1C end this facill- 
tatee its approximation to wmity where it appears. 


Considering the last term and non-dimensiondliting the derivetive .v lett= 


ag 3 We 
MS, 4= ag 7 > 


Oly pi s¥2 ae “st pe Bey one. ae 
16 2 dt ac) 2 dt V 
16-9 


Then ncen-dimensionalizing the entire tern, 


2 > ‘Be 
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2 sles ae Mig a rr” ( - 


Now ron-dimencionalizing the left side of the equation, after dividing bv 
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Thie equation may be rewritten in the more concise form, 


mt 
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There are then, the longitudinal equations of motion in nor-dimensicnal fort 


Se a ae “H i = 0 
M 73 M oy i Moi nal 7 
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These equations have been eeen to comoare term by tern, except for the 
inclusion of variation of thrust with velocity, with these developed in a 
more general sense. The acsumpticns and approximations have been cited. The 
dependent variables are u, @ and &. 

A more general case would be thet of assuming, instead of the homogeneous 
set above, the equating of the quantities above to forcing functione of the 
elevator displacement, §. This would require, in addition, en equation of 
motion involving moments about the elevator hinge line. This would be the 
case of stick-free stability and motion. 

It is noted that a striking difference between thege cquations and those 
more fam_liar to moet students is the form in which the etebility derivatives 
have resulted. ‘Since the problem is now ready for analysis of modes cf motion, 
damping factore and associated desired results, the queet’on arisee as to why 
this form of the derivativee and how they are to be evaluated. é 

fhe usual procedure is the separate evaluation of each derivative by 
Similarity to previous determinations or by wind tunnel or flight testing. It 
fs seen that certain of the derivatives in the three equations above are in 
combination. For a special but usual case of the solution of the equations of 
motion further combinations of these derivatives appear. A discussion of both 


the general case and this special one of the sclution foladowe. 





i) 


B. Solution of the Ecuations of Motion 

The most generel approach to the solution of the equations of moticn be- 
comes quite involved in that the solution is carried to the determination of 
an equation for each variable composed of both the complimentary function end 
particular integral. ‘See Reference 1. If the equations are restrictid to the 
stick-fixed case and the motion is aesumed resulting in the absence of en 
applied forcing function, the equations will be homogeneous and the solution 
resulting will be the complimentary function only. 

The discussion here will be further restricted in that the equations in 
terms of the variables will not be the end result This would merely give 
the particular motion resulting from certain initial conditions. 

The characteristic equation which results from the sclution of the dif- 
ferential equations being of the form, x = xe" , will be analyzed fori 
roots and information relative to inspection of these roots. This will show 
the character of the motion as regards stability, periodicity and demping. 

The nature of the characteristic equation, or etability quartic, will be 
set down from the determinant resulting from the solution, x = xe * ixX- 
perience has shown that the two sets of roote that form the sclution are 
usually of such widely separate magnitude that they may be separated. The 
Special case mentioned is based on this fact 

Since Scheudvel's solutions employing the quasi-static stability concepts 
are based on the premise that the roots are separable, the general characteris- 
tic equation 6111 ve set down without further comment and the special case 
will be streseed in preparation for the solution bv his quasi-static stability 
concepts. The errors involved in this approximation to the roots, especially 
as it increases with angle of attack, is discussed et length in Reference 1 


It is show that it ie a good approximation for lew angles of attack, l.e., ch. 
\ 


but becomes unaccectable for angles of ettack near ie especially for ngutre 1 
\ 


| 
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stability. 
l. The General Cage. > 
4 
Aesuming a solution of the form, x = x ars » where z isa réal or con- 


0 


ples constant, the indicated operations on the variablee are then carried out. 
The stability quartic résults from an expension of the following determinant 
that follows from the fact that the equations are coneistent if the determ nant 
Of their coefficients vanishes. 


Rearranging the equations by variables, 
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i > ns. | ne 
Sind lTodint | 
Be e+ Osi DS inv = : “- 
) 3 in C, | ok 
Sp ! “CL (235 in Vv yu bearc. sina | soe =e 
A. acs = MY a ae X-c =-c = er 
| 6 inv 2” ‘m2 Sn ae : 
i q i ce A i 
y . y y y : 
The Bcséumptions are made thet Cr gin 4s°- yv is negligible since near horizontal 


flight ic aesumed, and that apsroximately, To = De: 


oO 
The determinant of the cocffricients is then, using ae = sf ik 
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_ °°: Inv ~ 8 nv” el Bre 
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= C. (2% en j a5 
} EN ony oe “Log 
dC 
rm 2 & 2 2 2 2 
- k "eo. kk Le ae 
Tiny *y/ Sy SY a, By 7 omy By 
The quartic becomes, 
4 3 2 
ie. 1 4 Ls 4 + iF = 
| 1" .. a, a,- J. a 0 


Where , 
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2 m xy mM D‘d inv’ Oo inv Ly mY 


6 in Cy 
as “al Lote + sary | 


_9 = 2 
a Tox : 7 yas » o ln Cr 


ada Sl? ey! 


” - a in Cr C. 
C,(2 I OX 


— a 


oa , Cn C, (2 + See ) 
; (C + Ce ayes —— ~~ C||OFt— 
o,§ "i yy Lo 
Ad inc 
2 2 1G 
- C Ky | o2 (2 + ny 7 
GQ = ' 
ay eake Olea dC oi 
a _ - of (2 Ly 7 m Ke 
37 ‘sin Sin v 3B ft 
C(2 + sy 


The usual methods for the solution of quartics may be employed to solve 
this equatio as it stands. The constants have been arranged in the particu- 
ler fachion shown for reasons to be mentioned later. 

2. the Case of Distinct Sets of Roots. 

Experience has chown that the usual motion in flight consists of a long 
period, low frequency, lightly damped motion, called the phygoid or flight 
jasc: oscillaticn, and a short period, heavily damped motion, called the rotary 
oscillation. This suggests the factoring of the characteristic equation into 

two quadratic equations. 
The length of the period of the rotary oscillation hase been found to be 


at most of 1-5 seccnds while that of the phygoid usually ie of the order of 


30-60 seconds for most conventional planes. This suggests that, due to the 
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sil li duration of the rotary oscillation, the velocity is eesentially can- 
eth. Thue this mode is governed b~ equations A.5 and A.6 with the terms in 
u eliminated. The colution for the phygoid can neither ignore the change in 
meee ity or anvle of attack. FEowever, consideration of equation A.6 shows 
that, due to the lerge magnitude of the factor discussed in Appendix A, 

an approximation may be made that the terms containing this item are large 
compared to the other terms. Solving this equation for & , the terms in & 
in equatione 4.4 and A.5 are then expressed by a term in u. The phygoid will 
ve sOlved for, as regards the information mentioned previously, by considera- 
tion of the equations enclosed by the dashed lines in the equation array modi- 
fied by the approximation diecuesed above. The rotary oscillotion will be con- 
cerned by the cquations enclosed by the dotted lines in the array. 

et The Phugoid. 


Equation 4.G becomes, nelgecting terms not involving “ , 


oon a 2 
“Sinv “ya yA X= 0 
oY, 
OC 
Bem = (san emg 8 
Ecuations A.4 and A.5 in coefficient form are now, assuming a solution 
= ais 
©r the form, x = x .¢ F Ye 
m 
| oinD oi? 0 inv 7 
a+ 0, (so eee) ~ CO (———) + Ce 
mM P 
OC 
lea (in V) ee) 
> Yo ing Lx ¢ 
PX 
Which gives the quadratic, 
sys gt 
Bee 2 +a.Z+aq <2=0 
l Q 


where 











oo 
-o. (mb. gin) _, in) 
a, > *p‘S Inv’ @ inv Dey cm, 
and, 
oC 
a _=C (2 oon yy + C_C ( 2 =) 
fe) ) Fane ln V a6 c 
mR 
c#(2 + a “1 ye oe 
een in Vv - ' rm ie | 
- E “ay,  & InV in co 
mM oy ie (2 " = 
ae Oo inv 
The cusadratic has the roots, 
-a.t/a_ - ta 
Nt Ai fe) 
i a a 2 +" 
Zz 5 © Py gn it 1 
8) 
Where the damping factor, hy, is, By * the frequency or period factor is, 


i = Be 1 ae 
my Oo FF ; 


me ihe chort Period Oscilleticn. 
Neglecting the velocity terms, and equation A.4 in ite entirely, equations 


h.5 and A.6 become, in coefficient forn, 
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where, A, = (oc, -C. k ) and, A_= (e & ) k>.41.Cne root is, z= 0, 
iL ~H mo vi oO Mo sire] we ae K 
the others are, ¢ = R,, til, and thesdam@@ins fector is, BE, = - =, the frequency 
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facuor, 1. = ea FAY . 


‘Before any further evaluation of the quantities detcrmined uc to this 


Point, the quasi-static stability criteria are devoloped. 





C. CGuasi-Static Stability 

Before proceeding with the quasi-static stability, development, it may be 
helpful if a reiteraticn be made of the specialization aeeumed leading to 
thie development compared with methods of a more general nature ag in references 
land 4. Both cf these refcrencee discuss the epecieal cece of the approximation 
to the roots resulting from the assum tion thet the soluticn consists of two 
paire of complex rocte sufficicntly separeted in magnitude tc warrant a factor- 
ing of the charactericticse quartic. This facile solution is ueuil in most 
cages but may not be cound for weconvertional aircraft or even conventional . 
e@ircraft where additional facts must be consiiered. These are, flexing of the 
Bie cnent parte euch as wings anc tail eurfaces of the aircraft, whick would 
entail additional devrees of freedom, large changes in the coefficiente of the 
equations due to rapid depletion of a large fuel supplv and its effect on the 
derivatives. Many of these latter conciceratione are undergoing exhaustive 
Study at present and do not lend themselves to any concise generalizationa. 
From, say, a project engineer's point of view, especially in the initial etages 
of design, concern is given to information available from the most rapid end 
inexpensive data available. This is esvecially true et present eince epecifi- 
cations, the military in particular, require employmer.t of the ernie design 
throushout laryve ranges in altitude, speed and wei.ssht. The assumption of the 
presence of two quite widely separated modee, however, ie scund for mcst corven-: 
ticnal aircraft and for unconventional aircraft at certain phases of their 
flight history. The eonditionse that hold for the prugoid were seen to be a 
Slow orcillaticn of long period and weve lensth and relatively light damping. 


The rotary oscilietion ic a short, heavily damped oscillation. 


In regerd tc standard methods, the covfficients needed to determine the 
y 


constants for the final colutione cf tre equetions of moticn are obtained 


primarily from static wind tunnel tests. This is the rule for the values of 





ipl 


ecuetione 4.1 and A.2 ac seen in references 1, 3 and 4. The moment coeffici- 
ents are found by static etability consideretione for the variation of the 
nonent with & . See reference 5. A dynemic wind tunnel test is used to find 


C+. See references 1, h and 6. 
S| - 


The coefficients as they have been developed in thie paper will be evalu- 
eted in much the same manner as regards those concerned with linear forces. 
However, peeled of evaluation Of the logerithmic derivatives will be il- 
lustreted. As for the moment derivatives, the determination of will not de 
discucsed aseuning that it may be found by the usval methods. Tue steers longi: 
tudinal stability criterion, ae » is the quanitity that che next few pages wi. 
be concerned with. This quantity Dr. Schéubel determines in a quasi-static 
Way. Cauci-ctatic in that it is developed f.om assumptions made on the dynamic 
motion of the aircraft. The two concepts, quasi-static etability at equilibri- 
umn and et constant soveed, follow from the fact that an aircraft has two dis- 
tinct phases of motion following a disturbance. 

Cuasi-static stability at equilibrium of forces is essentially a considera: 
tion of the equilibrium conditions reached at a time in the flight history 
following a disturbarce at which the new velocity is obtained. This new veloc’: 
ty is the initial value plus the incrementel increase due to the perturoation. 
The cuaci-static stability at constant speed development holis for the rele- 
tively short period following the disturbance in which the velocity is assumed 
not to heave reaches ite final value. 

The dynamic response to an elevator deflection .s dieectec, so to speak, 
to g-ve these two quasi-static concents. The aircraft is initially in eteady, 
nearly horizontal straight flight. Ctarting from this steady etraight flight, 
we aftcuns that the elevetor angle, & » bac deen changed euvuddenly ba certa n 


amount, d§ , and we ack whet will hapven. 





ysis 


In tre iret tect, ae elevatcr displacement gives e« cheanse in longi- 
tudinal moment and a small change in lift. Thie Ictter chango ie so emll 
it is neglected. The change in moment dicturbing the equilibriun gives an 
ansular velocity abcut the lateral axis, and, due to this, the engle of at- 
tack is changed too. A chenge in the angle of attack givee rise to a chenge 
Grmerne Lift coefficient and by it @ change in lift. SO the equi libritm of 
forces perpendicular to the direction of flight ie €ieturbed too. All these 
effecte hold during the initial portion of tke period following the dietur- 
bance, i.e., during the first few seconds. In addition, it is assumed that 
the speed remains ocecntially coetent during this period. 

This then is the situation maintaining for quasi-etatic stability at con- 
etant speed and considers the change in force perpendiicular to the flight 
path caused by the initial curvature of the flight path and the change in mo- 
ments resulting from the dieburbance. 

Eventually the change in elevator ceusing the change in angle of attack 
thereby the lift coefficient results in the aircraft cttaining a new steady 
etate, that ie a new voint of equilibrium of forces and moments et a certain 
speec, V + dV, which is different from the initial one by cV. Cuasi-etetic 
stability at equilibrium cconcerne itegelf with thie portion of the flight hie- 
torv. From experience it is known that it takes an aircraft an appreciable 
time, ncermally several minutes, for the speed to adjust itself to a changed 


Clevator dicplacenment. 


iG 
\ 
1. Cvuasi-static Stability at Equilibriun. a 
\ Se 
The implication of an initial = an 
a LS 
Steady ctate of noticn necessitates oqui- a5 5 G7 . 


librium of forces and of moments. From 


Fizure 4 this ie eeen to be, 





Fisure 4. 
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Fos T-+D-W sin&= 0 


xX 
UF: L - W cos d= O 
1M: M=0 

IM, : mee 0 - 


Having started fxom nearly horizontal flight the flight path will remain 
nearly horizontal for sufficiently small changes so that the component of the 
weight perpendicular to the path remains almost unchanged. From the second 


equaticn, that of equilibrium of forces perpendicular to the path of flight, 





then, 
d(L - W cos & ) = dL = 0, since, ae before sin & = X= O. 
2 
v vo . av 
dL = dc —S+ 2c —S—=0 
~_ Ly 3 4) 2” Vv 
oc ac vc 
Bene: La 
SL Ox “oy “ame 
so that, 
Ss. ae ee 
sq aA oy (Saye 2) FI] PF 
or, 2 
oC, 
Sed inve - 8 — + dQ 
C (2 + Sa 
then finally, 3c, 
5.1 d inv = aX 
; dad x _ 0 inc 
: C.( 2am by 
L 0 inv 


wnere the subscript E indicates an equilibrium consideration with the agsump- 

tione implies above. This equetion is neither a pertial nor a total deriv«tive 

from a mathematical point of view and caution must be exercised in handling it. 
If the new state of motion is a steady one, which it must be, 

acy, = M- M, = 0. The moment coefficient depends on % , V, a and the angu- 

lar velocity, =< Since, however, it is assumed that the flight path hes no 

curvature at the new steady etate, and, indeed, in actuality euch would be the 


case, the angular velocity is zero. Thus the quantity, q, does not appear. 








There ie then, 
oC, 3c 3c 
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This then, is a new etability requirement thet applies to thet phase of 
the longitudinal motion described and implied from the previous discussion. 
This requirement certainly is sound. For an aircraft which does not have this 
quasi-etatic stability will be unstable in es much as for any disturbed state 
of motion, enfadrced by an elevator displacement, will show the tendency to 
move its elevator further in this direction, so increasing the deviation 


from ite initial state. 


dc 
, The usual static stability criterion, oe is considered either, 
ace dC dc dC 
~ oe LL or, - = - ) The quentit =" in either of the above is 
dC, dx"? ~ ax aa : ee? ee ear 


a& constant throughout the usual range of normal flight, i.e, wstalled flight. 
The value of x - x, can be found graphiclly from wind tunnel deta of Cy, vs. 
Cr as explained in reference 5. This is the more informative quantity since 
a shift in the center of gravity is more readily determined throughout the 
flight history. 

The firet term in equation C.2 is very similer to the static stability 
criterion above. However for comparison, equation C.2 is written in its en- 


tirety as, 
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where Caller is a new neutral point that takes into eccount variation of the 
neutral point with velocity, i.e., compressibility effects. Firet, it is 

noted that variation of lift coefficient with angle of attack differs between 
the two. This is due to the fact that the ordinary determination of the change 
of lift coefficient with angle of attack is accomplished in wind tunnel tests 
where the setreem velocity is kept constant. The results are then based on a 
constant dynamic pressure, whereas the precent development is not. The dif- 
ference between the two is about 1-30/o0 and is sumed negligible. The quantit~, 
as moves with respect to velocity, dynamic pressure, deflection of the fuse- 
lage, twiet of the wings and horizontal stabilizer. A comparison with the 
movement of the usual etatic stability neutral point will be made only with 
respect to variation with velocity and dynamic preseure. Figure 5 shows the 
variation with velocity of three different aspects of stability at sea level. 
These Be (e) static stability only, (b) quasi-stetic stability assuming the 


zero lift moment coefficient, C, , to be -0.02, and (c) quasi-static stabdili- 


M 
L 
ty with zero lift coefficient ° assumed zero. Figure 6 shows the comparison 


of the same quantities at altitude. Here the decrease in the stability merzgin 


(a) ®) (ay 
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for quasi-static stability is much less due to the presence of the density in 


the dynamic preseure in the denominator of the additional tern, i.e., Cy iy 


dy inV 
ad OX 


F ; 
2, Cuasi-Static Stability at Constant Speed. 

_ Recalling the discussion of the conditions holding following the imposing 
of a small elevator deflection on the steady state, it was seen that in the firs 
small interval following the change the speed does not change appreciably. Thue 
the quasi-static stability at equilibrium concept does not hold for this phase 
of the disturbed motion. It wae further seen that the equilibrium of forces 
perpendicular to the direction of flight is disturbed too. The equation of mo- 
tion for this direction shows what happens. 

1a PP ieice is exerted on a body moving in a certain direction, its path 
will ve oe See Figure 7. If an 
aircraft the motion may become rather 
complicated, for, due to the curvature 
of the flight path in a vertical plane, 
the component force of the weight per- 


pendicular to the path changes too. 





Restriction to small disturbances close pe 
to horizontal flight allows neglecting _ 
the changes in the weight component. : Figure 7. 


For instance, a change in angle of flight path of 15° givee only 3 o/o change. 


The centripetal acceleration is V a The ecuation of motion is therefore, 





W. ae : ; 
. a ap UL - W coed ) = al 
216 Oc 
aa, mee Bolla Ve 
ae =| 3 oe 3 28) | ¢ 2° 
at 
; ces cael oC, id te 
Tx Sq 8UP 5 








therefore, C ac 
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and, multiplying by —— and noting thet 4 = ples there is 
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Oc. iy 
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ice tern, » may be neglected since 50 ig small compared to Sar 


ee 


and, also, the usually lerge magnitude Of AA makes it negligible compared to 
unity. dhis is also the common approximation as seen by the neglecting of aa 


and Z in references 1, 3 and 4&. 
ac 
are Ll) 
Thus there is the reletion, Ae 2 Box oo 


tion of the forcee perpendicular to that of flight in the first inetant follow- 


» resulting from a coneidera- 


ae a distucbance. 

This curvature of the flight path or angular velocity has an influence on 
the equilibrium of moment about the lateral axis. So a consideration of tre 
Beiiitbriun of momente along this curved flight path involves the change in 
elevator deflection, the change in angle of attack and this angular velocity. 
The effect of velocity is omitted because of the aseumptions of one phase of 


the motion. There ie then, 





ac oy Ak , 2M ac aX 8H OF 0 
M ar dq dQ a ? 
or, 
OC, 
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This expreegsion is similar to the one found for the equilibrium of a steady 


State Of flight enforces by a certain elevator displacement. It contains tke 


fe $ 





stability tern, 


om —M m8 aa : 
ax [7 aie eq” do |, 


called "quaci-static stability at conetent speed”. A comparison of thie quan- 
tity with the usual etatic etebility is seen following a further evaluation. 
The angular velocity causes the angle of attack of the horizontal tail 


plain, X to change by the amount, 


? 
28 
ax = 2 but, ae * 49 75 
therefore, 


CO 7 


ig ae : 


This change in tail angle of attack gives rice to an increase in lift coef: 


My 





ficient of the horizontal tail plane, C. , and by it a decrease in moment about 


the lateral axis, a). 
dc. = + ac ag 
M Ls te 


and, having noted its dependence on Angular velocity, there ig, 


oc OC 
omen On a eta & (aby? 
dq 3X, aoe dq OX, hoo 


This quantity, Ou , mey be determined by a dymanic wind tumnel test de- 
q 


scribed in reference 1, or by use of the curved-flow technique referred to in 


reference 6. 
Oly dq 
The tern, ee , added to the etatic etability is always negative, 
q CK y 
Since Ou is a damping derivative. Thus this is an increase in etatic setabili- 
ty and the neutral point corresponding to this quasi-static stability ta be- 
ning aa one for ctatic etability by shout, 3) , Where, K= Cc 


this apparent improvement over wind timnel static stability is seen to be 


M Ses 


larger for-en aircraft of small wing loading at low altitude than for oe of 


the high wing loading et high altitude. 


The additional term added to the usual static et&bility neutre+ point 
looks eugpiciously like the ctick-fixed maneuver point, especially because of 
the setetement of the situation and the eseumptions made. This ig indeed the 
case and can be verified by comperiscn with the maneuver point foumd in ference 
L. : 

Concerniny the concepts of ieee static etadility, 1t is seen that both 
are imsrovenernts over Wind tunnel ststic stroility. However, three facte should 
be borne in mind concerning their effect end their use. First, both contain 
tme UEual stetic stability plus an acditional term ‘Secondly, their application 
Gemevce only hold fo. epecific conditions verculier to different phases of the 
Merent nietormyv- Inetly, for a specific aircreft they vary differently with Time 
pememetcrc, i.e., density, but all paremeters do not effect both. | 

Dr. Scheudsel pointed that in Germany during the lnst war certain of their 

tighters showod a marked improvement in longitudinal stability over that pre- 

dicted by wind? tunnel tects at mediun altutudes. Since these were of the 45C- 
500 mph class, thic was explainable by the abcve quasi-stetic etability coneid- 
erations. He 9@ia further, however, thet pilots complained of poor stadility 
@t high spoede at high eltutudes. In as much as doth concepts et ill predict en ; 
improvement, greater for equilibrium and lees for constent speed, the aeronseu- 
tree engineers were greatly puzzled. Further etudy acvoumted for this dy rea= 
son of flexing of the wings, the tail plane end the fuselage. The ac ountad! li- 
ty of theee phenomena in etability studies is very complex. For the lonzitudi- 
nal case a goce approximetion mace be obtained from grouping their overall 
effect ac a shifting of the neutral point. Tho ceriousness of their effect can 
be seen f.om mention of two simples agpects. F:ret, upon the sudden applicstior 
of an 6levetor defluction or a change in speed as in & gust, the elxsticity of 
present cecims at high speeds causes a twisting of the wing. Thie cruses the 


angle of attack tu vory epanwiee resulting in ao sh.fting of the center of 





pressure. Tho same applies to the horizontel tril eunaee. | secondly, @ 
shifting of neutral voint vortically f-oom flexing of the fuselage coupled with 
the very large drag forces at high epeeds proluces an appreciable destabilizing 
monent. 

That the problem ie complex analytically is easily seen from merely a 
superficial etudy of the present day litereture as, for instante, reference /. 
Dr. Scheudel pointed out the experimental difficulties in a description of an 
attempt in Germany to determine the deflection of a horizontal atabilizer with , 
speed. This involved a complex optical system employing an arrengenent of 
mirrors refliccting a slit light source the movement of which was recorded by a 
motion picture camera. In addition to this the test ee wae required to 


dive along a path determinsd by an ar.ay of floodlights to insure perfect lin- 


carity of flight path. 
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D. Fvaluetion cf the Perivativese. 
Sone of the quentities to be evaluzted are ijenticsl to those appearing 
in ctendard references. Mcthode for the determination of these will be men- 


tZoned but not discussed in deteil. 





Since all quantities appear in the general eclution, the items will be 


agiecuseed in the crder that they appear in the coefficients on perce 14. 


y cr . ectetic wind tunnel teete. 
A 
oy . kynemic winc tumnel teet by the whirling arm or oscillation mee 
q 


Gnoas mentioned in reference 1 or by the curved-flow technique referred to in 


& 
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reference G. 


a oo Ale letycor C~ ve. V can dve obtained from static wind tunnel tests 


ooen, 


DD, 

A Voc ae 1 
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Thre this vearics f..om a valuc of 2 


meee Mm US to about 3 fcr M= 1.0, 
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si Vy And ET ae . Theee aye ccn- yy, | 0 
Sidered toxscther because they are fcund . Se 
i es 
in muh the same manner. From winc tun- : ores al 
VioagM 
nel tests C. ve. V are fouw:.. iInygine Figure & 


manurectururs uveuzlly submit date of T ve. V fcr specific inetslletionse. These 
derivatives are then determined as in Fizures 9 end 10. It ie pointed cut that 


Figure 10 is plctted for a conetent anyle of attack. 
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Ficure 9. Figure 10. 
dc 


ss. A plot of Cn ve. V can be found from wind tunnel tests. The 


derivative is then found as shown In F’ gure 11. : _ 





Figure lic 





E. Evaluation of the Longitudinel Modes Using Cuasi-Static Concepts. 


1. The Phugoid. , 


The damping factor, k,, depends on the coefficient ay of equation B.2 and 





thie quantity is, OC, 
-c (2m _d8mt _¢ (o_in V) 
: “1° “p‘S inv Oo inv Dy °C, 
rons 
At relatively low velocity the approximate values of the derivatives are, 


aC 
3 2 - E, 2. 3 2 -1 (for a 4et) end Sea ie O. Thus the damping: is, 


one Q: S&S 3 a 
Ei, = 73 (¢,(3)| - = C, which is light. 


At higher speeds the effect of the last cuantity in R 


1 will be felt. This 


¢ 


will decrease the damping since it will overcome the increase due to the other 





quentities. 

The period factor, I; depends on doth the coefficients of equation 3.2. 
ime coetficiér.t, ay» using the quasi~etatic stability at equilibrium, equation 
usec, vecones, cf (2 . 0 In “Ly ac 

wera a 
O “Mae on.) 7 


The period factor nov becomes, 





Again, at low speede the approximations can be made that, 
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D2 X e 
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©o that, 
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3 in Cc, 3 Inc, 
5.2. dee Sa, 
oe Cr vV/2 Vv 
The period is then, - 
a. ot aoe v 
ee 0 ln a 6 
Vo *Sinv 
and the wave length is, 
ee) = 7 
3 ln Ce 
°*+ inv 
2a. The Short Period Oscillaticn. 
' ; el BY ae _ C 
The damping factor is, R, = - 2 Six Ca Kk). 
Approximate values for these quantities are, 
2 
= = -Q. sl ae ce. 
Cy Da, Oy OR and, M se 5 
% q 


ul 


The damping factor is the, RK, = 3 5 - (-5) = -5 which indicates heavy 


damoing. 


The ,»~criod or frequency factor was, 
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Fnd, using the values given before, it is seen that the item in the parenthesis 
is zerce, At leest the value in the parenthesis 1s very sml1l end it being to 
the second power, divided by four, in addition to the presence of the lorge 


eer, A ,» in the first term allows neglecting it in comparison with the firs 
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The pericd is then, T. = 7 = 





So then the period becomes, ; 





X ~ Xe Bp S) yo 


Fron the foregoing itis seen that depends on: : 
di 


meemeemter of gravity location: —===— . 
ee 5h 


Pen it ituie: 


Fe 
3( Wing loading: Most 
u) Size: hw 


5) Macs Concentration: 


5) Velocity: + 
Giveiccity: 7 
The decrement ise seen to be, 
ave 
(Cc -C, k) 
en BR, L M y 


0 = a 
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Although thie analysis is specialized and the modee of motion are of euch 








frequency and damping that they do not present much concern in design, the ap- 
plication of the quasi-static concepts has been chown. In regard to the im- 


portance of the modes, reference 8 pointe out that the importance of the phu- 


goid should not be m.nimized too much. 
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The evaluaticn of the coefficients of the non-factored quertic, equetion 
B.1, page 11, in terms of the queci-etatic et=bility concepts will not »e car- 
ried out. The coefficients have been arranged in 2 form fcr easy swvetitucion 
of theee quantities and a glance will show by whet degree they can be simpli- 
fied by such & eubdsetitrction. 

It is important, however, that a criterion for the separation of the quar- 
ieee into modés of slow ani fa@et frequency be establiehed. If thé ratic of (tae 
imaginery pert of the roots deser fin. ths phuzoid to thet of the short period 
Oseili@tion is emall compared to umity, then the aporoximetion is acceptedle: 
From the simplified frequency factors usec in the ont iow analysie, End, it 
i8 pointed out, the terms ret-.ined are the most important ones, it can be eeen 
that this critcriqn ic primarily dependent uvon the factors x - x, and As in 
the numerator of I. This follows since the other terms, althcuch they may vary 
eve rel#tively con&teant while the distencé vetveen the center of wravity end 


the neutral voint may be anywhere from 0.005 to 0.C5. The larve varietion of 


(¢4 


thre density factor is shown in Aopendix A. 
The oreccntation thus far has deelt with the etick-fixed cage only. The 


degree of freedom about the elevetor hinge line ic considered next. 
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FP. Introdwction cf the @levetcr Force. 

In Section A it ie recalled oe the equetion of lcngitutinel motion cue 
to the degrec of freddom abcut the elevator was locked, i.e., it inplied thet 
a chance of clevator angle was imorescet and heli feet. Now thie reetriction is 
removed and & rcevlt cbtsined in the form of a modification to the original 
Bolution of th@ ehcrt period oscillation. This is the stick-free case. It is 
interesting to notc here that the moet comncn texts on this subject treet this 
case quite differently. Feferencs 4 rulee crt env modification of the phugoid 
mode eni solvec the renaining equations, dropping terms in u from =n assumption 
of conetant speed, fo: two short perioc motes. Cne of the lstter is disregerdec 
@e the flapving cf the elevator about its hinge line, while the other is discus: 


1 


sed extensively ag the mode causing the "pornoising" motion. Reference 1 


yrites off the freeing of elsvators es a modification cf the modee ov © reduc- 
one, of the tail effect by reason of the wing downwesh and a possible contridu- 
tion to the velocity moment cerivative, M,, by reaeon of some constcnt moment 
due to elevator mass or 4 spring ettachment or friction. 

The equation of motion rosulting from this derzree of freedom consists of 


equating the moments resulting from freeing the elevatcr to the inertial tern, 


Oremoment cf inertia of the elevator times the enfulsar acceleration. 


° ArRCRAFT _- oo 
/ OG) ne eo 
L mpeman annem tin a + ~ aon ern ene retnie ey 6mm nerpnerere ee en emer een 0 en — 
ey = 
oa — 
cg HT eae hime 
™ @ 
Fopure ie. 
Thus, 0 
Fac tC 
i eel m i y of a M 
& “e 2 iS 





Css 
Ee) 


ae esp hs 3 ae ae on, Wee 7, S40 ED 


atm 

From the figure it is seen that a dow (positive) elevator deflection 
causee & nose dow (negative) noment about the aircraft center of gravity. It 
was neglected for the stick-fixed case but will now be included. Thies quantity 
i. in equation A.6. There are effects in directions consistent with the 
other equaticns also but these are agsumned small. 

It is also reasonable to assume thet there ise an increase in moment about 
the hintre line of the elevator due to the change in angle of atteck of the sein 
Age is pointed out in reference 4, this should be reduced to ae near poesible 
ae possible. A typical flghter of the lact World War had a value of Cho 25 Gye 
that of C, ng : aa may be made to be emeall by prover design and many methods 
of reducing it are discussed in references 4 and 9. The incremental angle of 
attack of the tail itself is emall for the portion of the flight histcry con- 
sidered here, and even that is coneiderably reduced bv wing downwash depending 
on the horizontal tail's position. The oe term will be neglected here, out 
certair.ly should be included if design prevents ite being made emall /fAny 


large value of C will detract from the pilot'e "feel" of the aircrgft. 


are 
Derivativee of the hinge moment with respect to retee of chanze of A, ond § 
aré nelsected also. : 

Ac in reference 4, any modificetion of the ghuzoid from feeeing the ele- 
vator will be ocsumed to be of no consequence. Thus the imprcvement to the 


étick-fixed approximation will only consider the short period oecillat on, 


which, as mentioned before, assumes constant velocity. Thus thie equation will 


& 
t 


neglect all changes with respect to velocity. 
The second term in Me is the coriollis force which reeulte from reference 

of the énguler motion about the aircraft center of gravity to the center cf gra 

vity of the elevator. The last term is the macs moment effect due to the angu- 


lar acclcration about the aircraft center of Gravity. 0 F.1 becomes, letting 
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then F.1]1 becomes, 
Pil t (lL A)F - t Babs t(1 + ok t SS 

Thic equation along with equations 4.5 end A.6 of Suction A gives tre eve- 
tem which shows the modes of motion for the stick-free caee. The complete s-re- 
tem of equations, including A.4, should be coneidered for a rigorove malig 
of the modes of motion. However, as was discussed previously, the echort peri- 


ed oscillation presumes conetant velocity and thie is again sesumed. Thue there 


is the following set of equations. 
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Ageuning a eclution of the forn, x = Xo e° , there ie the determine 1.t, 
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which is, regrouping, 
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or, eliminating the rcot, z = 0, and letting, f(z) ee z(C 


e ‘ 7 - 
and, f(z) e(l + a)(2 + ce ie abana is, 


- - 
F(z,t.) = f(z) te: f,(2) og aif ( 2) = 0 

If now the case is considered for the mase SEELCY » AA , large, which is 

usually the case except for gliders, then the quentit-, wee essentially vaniehe. 

So, approximately then, f,(z) = QO. The entire expression, F(z,t 2)? may be 


analyzed from the point of view of using the roots, 2,, of £,(2) in the evelue- 


ie 
ticn. 

The effect of freeing the elevators then will be the determination of the 
damping factor end frequency of the short period oscillation from 4 factor, 12, 
to be determined. 

Reviewing the eituation, the characteristic equation obtained from the de- 
terminant gives a fifth desree equatim in z. Elimination of the null root 
gives F(z,t,) which is fourth degree in z. Reference 4 then colves this equ- 


tion for two paire of roots which give the two short period oscilletions re 


previourly mentioned. 2 
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There is then, in mathemetical terms, 
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It is seen that, since both Shs and are negative, the real pert, which 
q 
indicatee damping, is negative and thus this motion is damped for this accep- 
tance of Sign. However, the value of Uy might go positive. For instance, re- 
8 


ference 9 ehows that a horizontal tail design using an NACA 66-009 airfoil eecti 
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with flap 0.30 chords in length, an air gap of C.CO5 chorde end with the hinge 


line at a positon of 0.50 flep chorde, reducee the U to zero betweer. 


Pox 
X = -9° ts & = ae while “hg is positive. Thie reference @lso ehows that 
C may go positive depending on the use or non-use of ahovnvalonce. The 


hs 


use of a hor balance, in general, appeere to mcke C s0 Slightly poelitive 


DAK 
&@lso, the magnitude veing very small 
Thus this mode of motion my be positively demped or possess moderste di- 


vergence. The frequency is given by, 
pep 
(1 2 ) sor | - Cy i 7. 
rs : 


2 f Toy ps ae 


where the term nezlected fe, (cy + a) , and ig congidered small oe in the 
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oe 
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previoue diecuesion. The terme in the breckete are negative, the last tern 
positive. The total value of the quentity is not too large. 

It is difficult to etrike e concrete canperison of the mode of motion 
eolved fcr here and thore found in reference 4. This solution does not eolve 
the equation by eae factorization but by the aesum tion of the quintity, tw 
being larre. In general, the rcots solved for by Dr. Scheubel are not the same 


as thoee caucing the corsoising motion described in reference 4, since there it 


is indicated that the nature of the porpoising mode is a function of Cy and 
oCL, a 
Shs - Here, it ic dependent upom Mae se and m, . It is not to be im- 


plies thet thie is not cuch a type of motion, rather thet it is not the sme 
im@plicitiv. 

Dr. Schéubel neglected the quentity ae , it fe dselieved, due to the fret 
that in Germany it woe the practice to design the vontrola euch treat it wea re- 
duced toa mnimun. He discussed the hondling qu: lities of > very large four- 


engined flying boet in which the stick forces per z were very emell without 


power boost, the exact valucs of which are not remembered. In regari to 
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with one foot on the rudd 
Inclusion of the C, term wor 
es wa 


- ¥ 
that it would detract or edd to both 
- ; “h , 
negative or pceitive in the retio «. fs Furthermore, in os much : 
wash would multiply the velue of the variable, & , which -ffects 


tor less than unity, this would further decrease the ratio itself. 
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G. Effect of Elevatr,r Impulse Tyoe Forcing Functiadn. 

Feiteratins, it is recalled thet the stick-fixéd cw#se involved the impor- 
tance of en €lovetor displecemert or a :vrust enc the ensuing motion we mnrly- 
zea in the lisht cf the elevet:r being locked in noeition. The etick-free cage 
allowed the elevitcr to floret freely abdovt ite hinge line. furthermore, the 
equations of motion were asfeumec homoseneous. An additional ccse ie now con- 
éidered In whicr the equations involve a2 forcing function end further this 
forcing Fimction is an elivator impulse, an elevator ceflection of short Gute. 
tion. This cace is analogous to the incorporation in the system of on on-off 
ty e autopilct. 

An elevator im ulse as show in Figure 13 ie imposed. Thie impulse ie 


giver. the cymbcl, Ig , and ig, 


or, 
ie t= J, Sat 


® 
a The first phace following the impulse, 


Gime c Of the short-period osciliation, is 


Conucidered and thue there are the equetione, Figure 13. 

; ho oc, 
A.5 x -5a A = Sa ee 
¥ Y - my +X -mK - = it oe 
A.6 x ne x gM Mg PAO nS A (7 ) 


Tre cffect of the im.ulee om the equztion of mction perpendicular tc the 
fina axis ic neglizible for conventional xircmft since the only force crerted 
is @ lift force, which, due to the short dummtia of the impulee, does not aot 
long erouwsh to beaffective. Thie ie not so fo taillees circreft, however, 
since for such aircraft the deflection of the elevotor, which ie »ctually a 
portion cof the main lifing foil trailing edge, changes the wing angle of Btteck 


and cauger an appreciable effect in thle direction. 





The effect of the impulse on the equition of moments causes on initial 
angular velocity which will be designated Z,, yut no angul-r displecement. 
If there ie no disturbance in flight path angle, which is nearly correct for 
conventional aircraft, this initial angular velocity will cmeist eseentially 
of solely a time rate of change of angle of attack. This is eo beczuse of the 
mass distribution of a conventional aircraft compared tc, sey, a tailless air- 


Ps ‘ N 
craft. Thus initielly, at © = 0; = F505 ASH 03 A= We 


Solving A.5 for §& , subetituting in A.6, there ie then for the short per- 


iod oscillation, 
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Imposin,; the first initial condition gives, A, = 0. Differentiating X then, 


1 
es 
: a i. 4a ~~ 
= = Lb 
x Woe (ya B, sin I,¢+ I, B, cos Ie) 
me 
The second initial condition gives, By = == . The Principle of Impulse and 
1 


Momentum stetee that the change in anguler momentum about cny axis equals the 


arpuler impuleéé of the applied forces about the axis, or, 
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Thig ¢valvatce the constant, 32 co tht consequently, 
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From equation A.5 thero is, xX = Cy ee ee sin I,¢. Inte- 
grating, 7 on 
a oe “a8 hy 
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if a ie eee 1 l 1 oe 
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This initial phace of the resoonee to an impulse tvpe forcing function is 
G6pictod qualitetively in Figure 14. In equation G.1 the vzlues of Ry and 1, 
ere the same ac thoce found in Secticn x” 

F., page 36. Thus the transient is 
heavily demned and typical values show 


it tc have reached tho cteady estate in 





oe-half periods. 


Examining further the eteady state 


Figure 14. 
of G.l, : 5 
yee, le Ong ty Pg Is 
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and letting, “thee ly = ce = ee there ig, oY = x If the de- 

nominator ic smell, the rec pmee is very seneitive to the impulse. Dr cSchen- 

pel citei an exanple of neutral etebility where the eteady state then becomes, 
C I : 

ms § 


eince, x= xX) 


4 eee. 
; an 
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low, pee anc, - C Using valuee of, 


2 
CY 
C2 
She 
&@ 
“y 
Ky 
as 
{I 
LS 
> 
‘a 
= i 
a 


9) 
UIT: 
fu 





4O 


30 
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which is appreciable for such a emil impulee. It is pointed out that thie 
“OK 

Ts 
teal horizontal tail area ratio of C.6 eccording to en empirically determined 


value of elevator cffcctivonesce, » correepon. s to en elevator area to to- 


curve in reference }. 


The seconc phase of the motion follcwing the impulse ie thot of the flight 
peth oecillaticn. Thut equetions A.4 and A.5 cre used. These should be modi- 


fied ae before for the variable, A . Fowever, it will be assumed that the de- 
oc 
rivative, sa ,» is approximetely zero and will be neglected. 


A.4 U- = ee ee ee a 
Sings 
A.5- C (2a Ju (f = fy (% ) 


Feroa the pect analysic of the two modes, the relL:tive mignitwies of their 
perioce, and the mawmitucese of their damping factore, it can de aesumed thet 
the regeponce te the impulse cf the rotary motion hee renched the etéady atate 
condition by the time the flight poth response tekes effect. Co the initial 

: n 7 — Y s = 7) 5 ra r— oO: 
G@@mcitions are, at, & = C: = @ , end us 0, an 
Again the fo.cin, function creatce no accountgble effect cither along the 


Wind axie¢ or perpendicular to it. Thus a Sood approximition ie to werk with 


honogencous equations. so, 
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The inclinetion of theflizht path, s* , may be defined, using Figure 15, 
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common time, © , eincc, t = ~ 
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This recult may be obtained by a slightly different Ipprocch’ Coneidering equr- 


Cr Ti Pee Soove, 











Differentiating C.2, 
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Ys -Y oe” = = ein I, 
2 
Therefore, 
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Uy ib 
u= he ‘ € sin I,¢ 
Coneldcring maximum u, 
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U Se ee 
max ee 
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and, neglecting enerfy loss due to de@mping, 
H + h* —= constant 
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eg 
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Which ls ae Sefore. 


Using the seme figurcs ac before, the m-xinmum guin in altityde very ap- 


proximately becomee, 
oO 


uh = $2Gtr 6.0995 = 95C feet 


eq 2 








Diecuesion end Vonclueions 


The quaci-eteatic ctabilitv concepts based on equilibrium of forces snd 
ments et varicus phacee of the flight history of en aircrsft follo.ing a die- 


turbance appear to be eomd for 2n acceptance of the description of the flight 


hietory steted sarlier. Neglecting adcitionel effecte such cs eloeticitr, t 


Gonceéyts account for the apvarent incre@se in etovility due ts compreseivilit 
éffectse associated with thé aircraft dr@s, lift an= momert If the voluet® of 
these quaci-etetic etabillty criterir couli be ~ecert ined individuilly and in 
their entiretv,  vresumebly they would inelude interference effects. As it aay 
no American literature imow to the writsr discussee lonsitudingl etedilitv 
compeariscns =t low ant high epyeede Dr. Schoudel C1 BCUS8ES the comperieon with 
the recultc chown ir F .wres 5 enf 6. This we epoarentl, b°sed m@m results of 
flight tests or the inclusion of the additional fectore in the qwei-etetic 
criterie anciytically dctermind. Thie letter woulé neccesarily neglect in- 
terv@erence effecte. 

For the special cace of seprration of the motion inte the phugoid and shor 
period modés xcod agreement with etendeard developments wes head neglecting the 
additional factore mentioned above. 


The mocificetion of the phugoil mode due tc Poeelnzg the elevators wee ne= 





elected ac ic ueveally dome. The refulting equrtion wee anc lyted for the sps- 
cial cace of the Bircvart meas deneity factor of ree menitude. This aesum- 
Peron wee uced to aboroximate only one peir of thé compley conftucate roots 
Comparieon with etandard references showed thie motion to be approxim=telv the 


Sporpo#Sing” mode althcugh not implicitly. 
Cf come intercet ie the Germen use of the quantity S-, square roct of the 


wing areca, ac the non-7imensiqalizing length term in convercion of forces and 
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moments to cocfficients. Actuollv, Dr. Schewrel did not point this ovt ce 
having wide Germen usace, but rather as being used ot Parmetedt Technical In- 
atitute. It m.ght be pointed out that the concept of cepect retio would pe 
lost by uee of thie lenswh. Thie might be better thin the nevuloue definition 
of mean aerodynemic cl.ord present in current Ancricon aerod-mamic literature 
For inctancc, refertnce 4, eupposeftly a fairly sccredited text, definee mesn 
eerodynamic chord as 1) an average chord, for a finite wing of varying section, 
Which, when multiplied by the average eection moment coefficient; dynemic pres- 
sure and wins area, gives the moment for the entire wing; 2) the chord of an 
imaginary airfoil which throughout the normil flight range hee the seme force 
vectors ag the three-dimenrional wing and cen be determined frem wind tunnel 
teste or a vzravhical analyeis which essumec c special shape to the wing plen, 
and 3) «a chord on which all thd forces and moments acting on thersurfece can 
be represented ec acting and on which the pitching mement coefficient 1# in- 
variant with lift cocfficient. Dr. Scheubel pointed out that the specicl ree- 
gon for the use of this quantity ws the fact that the ratio of the squcre root 
of the wing area to the tail length wae found to rance between ©.C end 1-25 and 
the avarase value found for some ninety aircraft we 1.01 The adventoge of 
approximating thie ratio to unity re@éulted in the edoption cof thie length. 

The cxtencion of euch @ method to the lateral or asymmetric case ie dis- 
Cuesed in gsencral terme in Appendix B. 

The development for the bngitudinal case ac a whole doee not consider roi 
the factore bearing on the problem by any means. It was pointed out by Dr 
ocheubel that thic was a good aporcximtion from the proicct encineer'e view- 
point and muet necessarily be improved as the deeign develope. For instance. 
tho effects of dowwash and the influence cf the powerplont have a lInrge ef- 
fect on the damping term. In any caec it is a elightly different view of the 


problem and ite worrhrests on that. 
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The mecs a a » as mentioned in reference 1 wis firet eug ested by 
Glauert for use in the non-limensionalizing of the equations of motion fhat 
it also follceowe from dimensional annalysia, e.g., the Pi Theorem, ie well kmown 
Due to ite extensive uec in the litereture, <« discu@sicn of it herw ie Only Wie: 
rantod by the slishtly different definition of it in this paper and the fact 
thet consideretione of ite magnitude for use in the approximtione in thepaper 
necessitate rccallinz wnat factors bear on ite magnitude. 

The difference in definition lies in the use of the mean wing chord es 
the represcentetive length in referencee 1 and 4. Furthermore, the mes or den- 
eity number uced in this paper is twice that used in theese references 

then is 
The gun tity year ae, 
2 wl? 4 ¥ 
An Me Be ‘5’ (=) 

from which it is seen to depend on altitude, Pp » on wing londing, W/S, end 
cnweight, YW. The tendencies with these perameters verying are as follows; 
a. Alitituce - As altitude yy increases. An altitude incresece from 
sea level to 40,000 feet caueces an increase in AA of frem 1 to 4 
b. Wing Loading - Since it is to the 3/2 power, an increaee of wing loading 
of 1:2 sivee an increace in A of Ii2neae 
c. Weight - Av weight increasce, ss decreases 2s (1/W)°. In thie cace a 1:2 
increase in weisht sive a 1:0.7 decrease in M . Typical values ga 

Lerge cliders at low altitude: » = 5-10 

Small fightere at high altitude: = 500 


Lerge commercial plenec: = 50-100 


APPINDIA B 


Discussicn of the Extension of Guesi-Ctatic Sti bility to the Lateral 


or Asynmet#ic Motion 


The possibility of the oxtensicn of the quasi-static stadility concepts 
LO the asynmhetric motion of ean aircra@ft seene to be remote insofar as tae 
same line of reasoning be followed. It will be recalled thet in both quasi- 
etaiic etability dcterminetions consideration wac made of forces in a single 
lincar direction and of moments about a point in one plane. In effect the 
ratio of changes in paramctere affecting ooth ceveteme lead to the quasi-sta- 
tic stability concepts. Cf fmdamental importance wac the fact that considér~ 
@tion of theee changes did not affect the additional linear force cquetton: 
“uch a correspondence does not hold in the latoral case.. 

Assuming a conventional aircraft it has been found that tne modes of mo- 
tion in the lateral case are four in number. Following the came line of rea- 
soning ag was done in the longitudinel case, & disturbance could, all facverm. 
ecisidered, cause a linear displacendént laterally, i.c., along the yu@xis, a 
rolling displacement, and a yawing displacement. References 4, 10 and 11 show 
wnat tne solution cf the lateral cquations results in a characteristic quair- 
tic which gives roots er Mop As, and hy, . In solving for these roots 
all neglect the presence of products of inertia, that is, they aesune that the 
orincipel longitudinal axie of the aircraft in ie line with the flight path, 


end stcume level flight or tan @ = 0. The moment composed of three components 
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Where the quantitice C ,» BP, «nd r are injicetive of sideelip, 2nguler veloci- 


ty cf roll, ond enguler vclocity cf yaw respectively. ihe coefficients Ff, 


1, ond 3, ere sependent «n initial concitione, on’y four veing rrsitrery, tly 


remainder are dctermined fiom trese. These indivicucl motions will be con- 
Sidered seperetely. - 

Comeidcring tho first ercup, the root, d,, ‘is emell end vewlly oositi 
giving light divergence. This is the su-called soiral motion cn? either e 
Gieturbance which ic pure sideelip, roll, or yew will develop on interdepen- 
dence of both rclling cnd yqawing momente on bcth sideslip angle end ycwing 


welocity. Follcwin; queci-etetic retecning, since the period is long, tie 





meticn indicatee t'e establishment of a unique rel tionship between initial 
disturbance incremer.te ené resulting forces and/or muments in equilibrium 
after a long period of timo. This does not follow cince either aileron, rui- 
ger or cidedglip, duc toa suet, Sav, Gey cance the moticn dnd 211 f3cece G7 
meversi moticn @re involved in the resulting motion. 

Tire eccond g@icup i@ characterized by ap Ionge root, A, peu lly neg: tive. 
Db, anc F. are emell in compezicon with P 


2 ra a 


and ie¢ cecentinlly pure roll. [To know reference corcider®s thie motion of anr 


eo tet the motion ie heavil- damped 


importerce. Reférencc 12 which includee the effect cf finite protuct of in= 
erevié, end inittel ircilination of the flight peth cxcludes «ny mention of 1 Beas 
Peferencc 15, which is an investigetion cf unectiefuctervy lsterel st vility 

Cn én sctuml curjcnt combat aircraft, does nct memtion it. ities evidence 
coupled with the fect thet the moticn involves only one decree of freedom 
ey docc nct lend itself to sny qmei-static stobilityv criterion. 


The third group is concemed with the all-important oscilletory motion 








ot 


* 


For fixed controls this coeiste of both a rolling motion and & yawing motion. 
It tco cen be induccé by ary ore or any conmoination of the possible lateral 
disturbances. 

rrom the foregoing it ic concluded thatumique representations in the quéesi- 
ebatic gence cannot ve deduced for the asyanetric motion. It appears from the 
comclucions cf reference 135 that awposeibvilitz gw ekt ve a censideration Of mae 
d@eree of fraom abovt the ruider hinge line. In the invéetigetion ovttined 
in this revort, ruider-fixed considerations were ruled cut in thet actual 


fligrt teets showed good aysreoment with predicted oscillations up to Mach num- 


Os 


ders correé ponding to a speed of 450 aph (density altitude 2t this speed was 


not cvecified). Insofar as a gradusl deterioration in @emying with speed 


recultins in ruddcr-fiee tects, a variation of the rudéer hinge moments, 


s 


beth Cc, _ and c ; 
he h5, 


this report that thie trend hes deen observed in high Mach number teste of 


3% 


with iech numoer is inferred. It is elco etated in 


cther contraol curfaces. Amplification on these intter mentionec tests is 
not upede. “wo facts must be pointec out here, however. First, the above men- 
tionoc. possibility is not in the quwmiegi-static cense of Dr. Schevvel's apvroach. 
Ccvecondly, reference 13 bascd its computed values at high speed on an extra- 
2Cletion of a cGetermination at an airspeed of 210 miles per hpur essumine e 
constant valuc of the etability derivetives throughcvt the speed range 

‘4 furthor conplication in the ectabliskment of quaci-etatic criteria fcr 
the lateral case sre the inclusion or exclusion cf the effects of products 
of inertie. Tneve effects ere diccueeod in detail in reference 1} and men- 
tioned in reference 12 for a specific aircraft confipurction. Thé effects 
of their inclusion was considered negligible in reference 13, but this cor- 


clusion muet be made weighed in the light cf the remarks mede cbdeve. 
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